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ABSTRACT 

We present an alternative technique for measuring the precursor masses of transient 
events in stars undergoing late stage stellar evolution. We use the well-established 
techniques of stellar population modeling to age-date the stars surrounding the site of 
the recent transient event in NGC 300 (NGC 300 OT2008-1). The surrounding stars 
must share a common turnoff mass with the transient, since almost all stars form in 
stellar clusters that remain physically associated for periods longer than the lifetime of 
the most massive stars. We find that the precursor of NGC 300 OT2008-1 is surrounded 
by stars that formed in a single burst between 8-13 Myr ago, with 70% confidence. The 
transient was therefore likely to be due to a progenitor whose mass falls between the 
main sequence turnoff mass (12 M Q ) and the maximum stellar mass (25 Mq) found 
for isochrones bounding this age range. We characterize the general applicability of 
this technique in identifying precursor masses of historic and future transients and 
supernovae (SNe), noting that it requires neither precursor imaging nor sub-arcsecond 
accuracy in the position of the transient. It is also based on the well-understood physics 
of the main sequence, and thus may be a more reliable source of precursor masses than 
fitting evolutionary tracks to precursor magnitudes. We speculate that if the progenitor 
mass is close to 17 Mq, there may be a connection between optical transients such 
as NGC 300 OT2008-1 and the missing type II-P SNe, known as the "red supergiant 
problem." 

Subject headings: galaxies: individual (NGC 300) — galaxies: stellar content — stars: 
evolution 
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Introduction 



Energetic and luminous transients mark the final stages of massive stars' evolution. These 
transients include core-collapse explosions leading to Type II and Ib/c supernovae (SNe), luminous 
blue variable (LBV) outbursts, and other more myst erious transients su ch as the optical transients 
in NGC 300 (herea fter NGC 300 OT2008-1, following liter ger et alJlioOfll ) and NGC 6946 (SN 2008S; 
Prieto et al.ll2008al ). Current stellar evolution models indicate that the mapping between progenitor 
stars and the type of supernova is determined primarily by the zero age main-sequence mass and 
the star's mass-loss history. The theoretical expectations are that the most massive stars (> 25 
Mq) experience tremendous mass loss at the end of their lives, expelling their hydrogen and/or 
helium envelopes before exploding as SN Type Ib/c; somewhat less massive stars (~ 10-25 Mq) are 
assume d to lose less mass a s they evolve, and e xplode as SN Type II wi th their hydrogen envelopes 



intact (IWooslev et al 



2002 



Heg er et al J 120031: iLim ongi & Chieffi 200 71). This picture is supported 



by ob servations ([ James &: Andersonll2006l : iPrieto et al.ll2008bl : lKellv et alj|2008l : I Anderson &: James 
2008). On the other hand, theoretical expectations for the progenitors of optical transients such 
as N GC 300 OT2008-1 and SN 200 8S are lacking. Empirically, LBVs are expected to be > 20 
M dSmith et alJbooi ISmithl hoOTh and estimates for NGC 300 OT2008 - 1 and SN 2008S place 



the progenitor mass somewhere between 6 and 20 (jPrieto et al. 



2008a 



Thompson et al 



2008 



Smith et al.1 bpOfll : bond et al.1 bood : botticella etalThoOflh . However, it is unclear whether the 



latter are SNe resulting from the collapse of Q /Ne cores (hereafter refer red to as "electron-capture 



SNe" ; IPrieto et al 



2008a 



Thompson et al 



2008 



Botticel la et al.ll2009i ). an extension of the LBV 



outburst phenomena to lower masses (jSmith et al.l l2009l ) , or a new class of transient due to an 
unknown mechanism. 

At this juncture, progress in understanding the origins of these transients would be aided by 
the addition of observational constraints. The increasing prevalence of time-domain observations is 
adding to the number of transients identified, but the most basic property of the transients — their 
stellar mass — remains largely unknown. 

One appr oach to constrainin g the masses of the stars showing transients is to use precursor 
imaging (e.g., ISmartt et al.ll2009i ). If images of a star exist before outburst, stellar evolutionary 
tracks that pass through the color and magnitude of the likely precursor can then be used to 
estimate the star's initial main-sequence mass. Core-collapse explosions which disrupt the star 
require precursor imaging to identify the progenitor, whereas imaging after the transient is possible 
for LBVs and other transients that do not completely disrupt the st ar. This second class of objects 
can be misidentified as SNe and has been labeled "SN impostors" (jVan Dvkl 120071 ). 



Several groups have noted two unusual transients discovered in 2008 with spectra sim i lar to 



Type Iln SNe, but whose luminosities are much l ower than typical S Ne (jStanishev et al 



Bond et al. 



2008 



Monard 



2008 



Bond et al 



2009 



Berger et al 



2008 



2009). These two transients- 



SN 2008S and NGC 300 OT2008-1 — have no optical precursor, but do show a very red source in 
the mid-IR bands covered by IRAC and MIPS on the Spitzer Space Telescope (jBerger fe Soderberg 



2008al lbl: IPrietol bood : IPrieto et al.1 l2008al : bond et al.1 bood : feerger et all bood ). These two pre- 
cursors have been interpreted as defining a new class of heavily dust enshrouded s uper asymp- 
totic g iant branch (AGB) stars, which may have exploded as ele ctron-capture SNe (jPrieto et al 



2008al : iThompson et al.1120081 : Botticella et al.ll2009h . Alternatively, ISmith et al.1 (|2009j) suggest that 
SN 2008S is not a core-collapse supernova, electron-capture or otherwise. They instead propose 
that SN 2008S is a "SN impostor" with a super-Eddington wind emanating from a star with mass 
> 15 M & . 

Recent light curve and spectroscopic observations of N GC 300 OT2008-1 confirm its overall 
similarity to SN 2008S ( Bond et al.ll2009l ; iBerger et al.ll2009l ). The spectra show evidence for both 
a wind and fallback of material from a previous mass-loss event. However, due to the uncertainties 
in stellar evolution models, translating bolometric luminosities to mass estimates results in a fairly 
broad range of possible masses for these stars. 

While direct imaging of precursors clearly has promise for estimating masses, it suffers from a 
number of substantial limitations. First among these is the requirement that the precursor imaging 
actually exists. The majority of past transients have neither pre-existing Hubble Space Telescope 
(HST) imaging, nor sufficiently accurate astrometry. The second major limitation is that even 
when precursor imaging is avail able, interpretation of that imaging depends on the most uncertain 
stages of stellar evolution (e.g., iGallart et al.ll2005l ). Existing studies typically estimate the mass 
of a precursor by fitting evolutionary tracks to its color and magnitude. However, uncertainties in 
stellar evolution models related to rotation, mass loss, pulsation, internal mixing, the formation 
of dust in stellar winds, and convective instabilities in shell-burning layers all contribute some 
difficulty in determining the mass of a single highly-evolved star (see discussion in £ j4,3.5p . 

In this paper, we demonstrate an alternative approach using stellar populations. Our method 
is based on the fact that most transient events occur in the last stages of stellar evolution and 
within groups of stars that share a common age and metallicity. Thus, even after an individual star 
has evolved, the remaining stars still provide information about the age of the stellar population 
that hosted the transient, even if the source of the transient is no longer visible. Once that age is 
known, one can infer the mass of the star that exhibited the transient, since it likely corresponds 
to a star that has recently turned off the main sequence. 

This meth od has been used to age-date SNe for which the surrounding stars could be resolved. 
Efremovi (jl99ll ) identified a young cluster surroundin g SN 1987A and est imated its age as ~ 10 
Myr; this estimate w as later refined to 12 ± 4 Myr (jWalborn et al.l Il993l ) and then 12 ± 2 Myr 
(jPanagia et al.ll2000l ). IVan Dyk et al.l (jl999l ) studied stellar populatio ns in the vicinity of a numbe r 
of S Ne and were able to estimate progenitor masses in a few cases. iMaiz-Apellaniz et al.l (|2004) 
and IWang et al.l (|2005l ) estimated the age and mass of SN 2004dj using not resolved stars in its 
originating cluster but by finding the combination of spectral energy distributions (SEDs) which 
best fit the integrated cluster light. They arrived at age estimates of 13.6 Myr an d 20 M yr, 
respectively, which correspond to progenitor masses of 15 and 12 M@. IVinko et al.l (|2009l ) fit 
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isochrones to the color-magnitude diagram (CMD) of resolved stars in this cluster and placed the 



age at ~ 10-16 Myr, corresponding to a progenitor mass of ~ 12-20 M & . ICrockett et al.l (j2008) 
attempted to date SN 2007gr using SED fitting, but found two possible solutions at 7 ± 0.5 Myr 
and 20-30 Myr. 

Our method takes advantage of the fact that most stars form in stellar clusters with a common 
age (At < 3-5 Myr) and metallicity. Indeed , over 90% of stars form in rich clusters containing 
more than 100 members with M > 50 M Q (|Lada &: Ladal 120031 ) . The stars which formed in a 
common event remain spatially correlated on physical scales up to ~ 100 pc during the 100 My r 



lifetimes of 4 M & stars, even if the cluster is not g ravitationally bound (IBastian Goodwinll2006l ); 
we have confirmed this expectation empirically in lGogarten et al.l ( 20091 ). Thus, it is reasonable to 



assume that most young stars within approximately a hundred parsecs of a massive-star transient 
are coeval. The age of a transient's host stellar population can then be recovered from the CMD 
of the surrounding stars, using well-established methods for deriving the star formation history 
(SFH). 

Herein, we analyze the specific case of NGC 300 OT2008-1, by extracting stars from a small 
aperture around the locat ion of the transient, u sing photometry from the ACS Nearby Galaxy 



Surve y Tre asury (ANGST; Dalcanton et al.ll2009h . Using the methods described in I Williams et al. 
( 2009 ) and Gogarten et al. ( 20091 ). we solve for the SFH in the region of the transient. We find 



a very well-defined burst of star formation that occurred 8-13 Myr ago. The burst allows us to 
accurately age-date the evolving stellar population which produced the transient. The mass of the 
transient can then be linked to the masses of stars that have recently turned off the main sequence 
but have not yet exploded, which in this case is 12-25 M@. We will show that with this method, 
it is possible to establish firm limits on the mass of the precursor, using only imaging taken after 
the transient itself. 

The outline of the paper is as follows. In $2] we discuss the stellar photometry used to analyze 
the age of the transient. In §3.11 we discuss the recovery of the SFH. In §3.21 we calculate the 
resulting stellar mass of the main sequence star which evolved to produce the transient. In £@]we 
discuss the implications of our results, strategies for expanding the use of this method, and possible 
limitations. We conclude with £j5j 



2. Data and Photometry 



The location of NGC 300 OT2008-1 (a = 00 h 54 m 34.552 s , S = -37°38'31.79" (J2000 UBerger Soderbere 



2008bl ) is contained both in one of our ANGST fields observed 2006 November 8, and in an archival 



observation taken 2002 December 25. The ANGST data had exposure times of 1488s in F475W, 
1515s in F606W, and 1542s in F8UW. The archival data had exposure times of 1080s in F435W 
and -F555H^, and 1440s in F814W. Each of these observations was split into two exposures, which 
were calibrated and flat-fielded using the standard HST pipeline. Photometry was done using 
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DOLPHOT, a modified version of HSTphot (lDolphinll2000h optimized for ACS. Det ails of the pho 



tomet ry an d quality cuts used for t he ANGST sample and archival data are given in I Williams et al 



(|2009l ) and iDalcanton et al.1 (|2009l ). We use only the highest-quality photometry for this analysis. 
Although the ANGST data was in fact taken befo re the transient occurred, there is no star visible 
at the location of the transient (|Bond et al.ll20 09). and thus the precursor was not visible at optical 
wavelengths. 

Artificial stars, used to characterize the completeness of the photometry, were inserted and 
detected using DOLPHOT. We have approximately 200,000 artificial stars for each field in the 
region in which both the archival and ANGST fields overlap. We use fake stars only for the overlap 
region so the effects of crowding will be identical for the two fields. 

We select stars in a 5" radius around the transient location (Figure H]). This translates to 
~ 50 pc at the distance of NGC 300 (2.0 Mpc), which is a compromise between including as many 
coeval stars as possible while limiting contamination (see M4.3.4I for a discussion of the choice of 
selection radius). The CMD shows a truncated main sequence, indicating that star formation in 
this region terminated more than ~ 5 Myr ago. However, the CMD shows no core helium burning 
stars, suggesting that there was no significant star formation during the interval between 25 and 
300 Myr ago. Thus, even without a detailed analysis, the CMD alone indicates that the transient 
was due to a massive star that formed during a burst of star formation between ~ 5 and 25 Myr 
ago. In the next section, we analyze the SFH in depth, giving tighter constraints on the age. 



3. Analysis 
3.1. Star Formation History 



Overlaying isochrones on a CMD is the traditional method of determining the age of a young 
stellar population; however, this method relies on only a few main-sequence stars to establish age. 
In contrast, fitting the entire CMD uses the full wealth of information and can give tighter age con- 
straints. Comparing the observed CM D to a set of model CMDs to derive the star formation history 



(SFH) is a well-established technique ([Hernandez et al 



19991 : iDolphin 



2002 



Skillman et al 



1999 



20031 : lHarris fc Zaritskv 



Gallart et al 



2004 



1999 



Gallart et al 



Holtzman et al 



200.51 ). While there 



are many different codes available, the basic procedure is the same for all: stellar evolution models 
are used to predict the properties of stars of different masses for a range of ages and metallicities. 
From the predicted luminosity and temperature, the magnitudes of the stars are determined for a 
given filter set. For each age and metallicity, stars are placed on a synthetic CMD in proportion 
to an assumed initial mass function (IMF). These CMDs are then linearly combined, with distance 
and extinction either fixed or included as additional free parameters, until the best fit to the ob- 
served CMD is found. The ages and metallicities of the CMDs that went into the best fit tell us the 
ages and metallicities of the underlying stellar population, while the weights given to the CMDs 
provide the SFR at each age. 



-10 12 3 

F606W-F814W 



Fig. 1 — Left: subsection of the 3-color (FA75W, F606W, F8UW) 2006 ACS image of NGC 300. 
The location of the transient is marked with a star, while the circle indicates the 5" radius from 
which we selected stars. North is up and east is to the left. Right: CMD of stars in the selected 
region. Photometric errors as a function of F814W magnitude are shown on the right. Isochrones 
corresponding to th e edges of our age bins for SFH derivation are plotted. All isochrones are from 
Marieo et aD (fcoosh and have [M/H] = -0.4, Ay = 0.1, and m - M = 26.5. 
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To derive the SFH for NGC 300 OT2008-1, we use MATCH, which finds the maximum- 



likelihood fit to the CMD (|Dolphinll200a ) . We assume an IMF with a slope of -2.35 (|Salpeterlll955l ) 
and a binary fraction of 0.35. MATCH only allows a single value for the slope of the IMF, but 
given that our CMD only includes stars with masses > 1 M@, adopting a single Salpeter slope is 
likely to be a valid a ssumption. Synthetic CMDs are constructed from the theoretical isochrones of 
Marigo et al.1 ( 2008) for ages in the ran ge 4 Myr-14 Gyr. The isochrones younger tha n ~6x 10 7 yr 
were a dopted from Bertelli et al. ( 1994 ). with transformations to the ACS system from lGirardi et all 
(|2008l ). Age bins are spaced logarithmically since the CMD changes much more rapidly at young 
ages than at old ages. The edges of the age bins shown in this paper are 4 Myr (the age of the 
youngest isochrone), 8, 13, 25, 40, 63, and 100 Myr. Additional age bins go back to 14 Gyr. 

Metallicity is constrained to increase with time within the range —2.3 < [M/H] < 0.1. As 
additional free parameters, the distance modulus is allowed to vary in the range 26.3 < m — M < 
26.7 and extinction can vary in the range 0.05 < Ay < 0. 50. Also, up to 0.5 mag of differential 
extinction may be applied to young stars (< 100 Myr). The lSchlegel et al.1 (|1998l ) value for Galactic 
extinction is Ay = 0.042 in the line of sight to NGC 300, but we expect the total value to be higher 
due to local extinction within NGC 300 itself. 

Comparisons to model CMDs were carried out within bins of width 0.1 mag in color and 0.2 
mag in magnitude. For a small number of stars, choosing bins that are too small results in so few 
stars in each bin that the accuracy of the fitting suffers. Our choice of bin size reduces this problem 
while ensuring that the number of bins in the CMD is substantially larger than the number of free 
parameters in the fit. 

Completeness and photometric biases were accounted for by including the results of the artifi- 
cial star tests. We supplied MATCH with the input and output magnitudes of the artificial stars and 
whether they were detected above the quality cuts of our photometry. We include only the portion 
of the -F606iy + F&1AW CMDs that were more than 50% complete, which corresponds to magni- 
tudes of F606W < 27.9, F8UW < 27.0 for the ANGST data and F555W < 27.6, F8UW < 27.0 
for the archival data. We also recovered the SFH using the F435W/F475W + F814W filter combi- 
nations and found consistent results, but since the depth in the B equivalent filters is less than in 
the V equivalent filters (^5551^/^6061^), the error bars on the resulting SFHs were considerably 
larger. 

The resulting SFH of this region over the past 100 Myr is shown in Figure [2j The SFH shows a 
single star formation event in the range 8-13 Myr ago, with no significant additional star formation 
in the age range of 0-100 Myr. This isolated burst is the likely source of the progenitor. The time 
of the star formation event agreed between the ANGST-derived and archival-derived SFHs, with 
the amplitudes of the star formation events agreeing within the error bars. We note that many of 
the fainter stars in this region were formed at a low SFR over the past 14 Gyr; however, we plot 
only the past 100 Myr to highlight the event that likely produced the progenitor. 



Error bars in Figure [2] are the quadrature sum of the errors from uncertainties in distance 
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and extinction and the 68% confidence interval from Monte Carlo tests, which assess uncertainties 
due to limited numbers of stars. Monte Carlo simulations were run by sampling stars from the 
best fitting model CMD determined by MATCH. These stars were then given as the input to 
MATCH, and the resulting SFH was compared to the SFH from which they were drawn. This 
process was repeated 100 times, and the scatter in difference between the input and output SFHs 
was incorporated into the error bars in our reported SFHs for each region. These trials indicated 
that the residual recent star formation at ages older than the burst was unlikely to be real, and 
that the SFH within the past 50 Myr could be explained entirely with a single burst. Monte Carlo 
simulations assess uncertainties due to Poisson sampling of underpopulated regions in the CMD, 
but they are not sensitive to systematic uncertainties in the models themselves. Hence, our error 
bars for the SFH do not include model uncertainties, but should be viewed as the range of possible 
values for this particular set of models. We discuss model uncertainties further in section §4.3.51 



Extinction was found to be Apeoew 
for NGC 300 of E{B - V) 



0.09 ± 0.04. 



Gieren et al.1 (120051 ) find a total reddening 
0.096 ± 0.006 mag, w hich corresponds to an extinction of Ay = 0.3 



assuming the extinction law of ICardelli et al.1 (|1989i ). Given that the extinction in NGC 300 varies 



with po sition (Roussel et al. 20051 ). it is not surprising that the value for this particular region is 



smaller. iBerger et al.l (|2009i ) estimate the extinction of the transient itself from SED fitting (after 
much of the dust surrounding the progenitor may have been destroyed) and find Ay = with a la 
upper limit of Ay < 0.6. Our value is consistent with this estimate also, although we note that the 
value we derive is the average extinction of the stars around the transient and not of the transient 
itself, which would i nclude the ci r cums tellar dust that prevented the detection of the progenitor at 
optical wavelengths. iBond et al.l (J2009J) find from their spectroscopy that the transient itself could 
be reddened by as much as E(B — V) ~ 0.4, corresponding to an extinction of Ay ~ 1.2. 

The best-fit distance modulus was found to be 26.41 ± 0.09, which a grees with the yalue o f 



26.50 found using the tip of the red gia nt branch for the entire ACS field (IDalcanton et al 



as well as previous distance estimates (I Butler et al.l |2004| ; iGieren et al 



2005; 



Rizzi et al 



2009), 



20061 ): 



however, this is by construction given the limits placed on the distance modulus. 

The current metallicity returned by MATCH is [M/H] = —0.44 ± 0.17, which is in agreement 
with the value of [M/ H] = —0.27 ± 0-1 5 . at th is galacto centric radius (2.7 kpc) from the observed 
metallicity gradient of iKudritzki et al.l ([2008). To test the effect of metallicity on our results, we 
fixed the metallicity at [M/H] = —0.7 ±0.1 when deriving the SFH. The SFH still shows a burst at 
8-13 Myr ago, though with an increased probability of the burst beginning earlier (20 Myr ago), 
presumably due to redder stars being interpreted as older rather than more metal-rich. Forcing the 
stars to be at solar metallicity results in a burst that is too young, 4-10 Myr ago, but that scenario 
is clearly inconsistent with the measured metallicity at this radius. 
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3.2. Mass of Progenitor 

The relationshi p between stellar age and the main-sequence turnoff mass is plotted in Figure [2] 



(jMarigo et al.ll2008l ). From this plot we can infer the progenitor mass for a stellar population of a 
given age. The star formation event which gave rise to the transient happened between 8 and 13 
Myr ago, which corresponds to a main sequence turnoff mass that falls in the mass range of 12-17 
Mq. Assuming that the transient was due to an evolving star, the transient precursor had a mass 
higher than this lower limit. 

To determine the uncertainty in our turnoff mass estimates, we ran a series of Monte Carlo 
tests using artificial star formation bursts of varying amplitudes and durations and measured how 
well these SFHs could be recovered. We varied the amplitude of the bursts (i.e. SFR or number of 
stars formed in a constant time period), the start and end times of the burst (for varying duration), 
and the age of the burst for similar durations. Each simulation was repeated 50 times. 

The amplitude trials were set so as to produce a certain number of upper main sequence 
stars (defined as -F606PF — ^8141^ < 0.3 and FQ06W < 26). Trials were run with numbers of 
approximately 25, 50, 100, and 150 (exact numbers varied slightly with each trial, since CMDs were 
created by sampling from the fixed SFH). The burst duration was set to 8-13 Myr to match the 
derived SFH of NGC 300 OT2008-1. Results are shown in Figure [3l If each burst were perfectly 
recovered by MATCH, the cumulative distribution would go from to 1 within the gray shaded 
region which shows the burst duration. All of the trials have a tail to older ages and smaller masses, 
indicating the tendency of MATCH to interpret stars as older than they really are, but this accounts 
for less than 40% of the total star formation in all cases with 50 or more upper main-sequence stars. 
The star burst with ~ 50 upper main-sequence stars is the closest match to the derived SFH, which 
had 55 upper main-sequence stars. We see that the burst is not well constrained with only 25 upper 
main sequence stars, but it is at 50 and above. We conclude that at least 50 upper main-sequence 
stars are required to use this method. 

We next set all duration trials to produce approximately 50 upper main-sequence stars, and 
then varied the starting and ending times of the burst. The end time trials began at 13 Myr ago 
and ended at 5, 6, and 10 Myr, and the start time trials began at 10, 16, and 20 Myr ago and 
ended at 8 Myr (Figure 0]). These can be compared with the 8-13 Myr burst from Figure [3l We 
see that very short duration bursts (2-3 Myr) are not able to be recovered, with the resulting SFHs 
extending to significantly younger and older ages. Longer bursts are better constrained, especially 
at the upper turnoff mass limit. Despite the long tail to older ages, this actually represents a small 
range in absolute turnoff masses, as the 7-10 Mq range stretches over almost 20 Myr. 

Since we used logarithmic time binning, we could not fix a constant linear time duration for a 
variety of ages, so instead we chose the time bins which gave durations as close as possible to the 
derived duration of 4.6 Myr. Figure [5] shows the results of these trials, with start times of 4, 5, 6, 
10, 16, and 20 Myr ago. From these plots we conclude that bursts of this length are more easily 
recovered at more recent times. Bursts that ended more than 10 Myr ago are frequently interpreted 
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as ending more recently. We speculate that this problem may be related to poor sampling. If upper 
main-sequence stars are present, the CMD must be young, but the absence of these stars could be 
due to chance at the low levels of star formation that we are probing. Bursts with very recent end 
times show no tails to more recent ages because the youngest isochrones in our models are at 4 
Myr. 

We use a subset of these tests to assess the uncertainty in our mass estimate. The question we 
are asking is whether a burst that started earlier or ended later might be interpreted by MATCH 
as the 8-13 Myr burst we derived, and with what probability this might occur. For the start time 
variation, we define a "matching" burst as one that has a SFR < 1(T 5 M yr^ 1 for ages older than 
13 Myr and a SFR > 10~ 4 Mq yr _1 at some point in the 8-13 Myr range. These cutoff values 
are taken from the derived SFH of NGC 300 OT2008-1, which has SFR > 10 -4 M Q yr" 1 during 
the burst and < 10~ 5 Mq yr _1 at other times. For the end time variation, we likewise define a 
"matching" burst as one that has a SFR < 10~ 5 M Q yr" 1 for ages younger than 8 Myr and a SFR 
> 10" 4 M Q yr" 1 at some point in the 8-13 Myr range. We ran 50 additional Monte Carlo tests 
for these bursts, for a total of 100, but the results remained the same as for the original 50 runs 
testing the recovery of burst duration. 

Of the bursts that started earlier (8-16 Myr and 8-20 Myr ago, top center and right panels 
in Figure H]), we found that the probabilities of recovering a single 8-13 Myr burst were 40% and 
22%, respectively. This means that if we derive an 8-13 Myr burst for an observed CMD, we can 
be 60% sure that the burst started after 13 Myr [M prog > 12.4 Mq) and 78% sure that it started 
after 16 Myr (M prog > 11.9 Mq). Of the bursts that ended later (6-13 Myr and 5-13 Myr ago, 
bottom left and center panels in Figured]), we found probabilities of 30% and 26%, respectively. 
Thus, we can be 70% certain that the burst ended by 8 Myr ago (M prog < 17 Mq) and 74% certain 
that it ended by 6 Myr ago [M prog < 20 Mq). Based on these tests, we report our uncertainty on 
the 12-17 Mq turnoff mass estimate as 70% (~ la in a Gaussian distribution), setting the lower 
limit to the transient precursor mass. 

At the higher mass end, the mass of the precursor is limited by the most massive post-main 
sequence star that is expected to exist for a burst of star formation between 8-13 Myr ago. As we 
show in Figure [61 isochrones for 8-13 Myr old bursts show no stars currently more massive than 
~16-25 Mq, with ~1 Mq variation between different stellar models. If these models are correct, 
then the progenitor mass lies above 12-17 Mq but below 16-25 Mq. 

To test the effect of varying the IMF and binary fraction on our results, we generated artificial 
SFHs with a burst from 8-13 Myr, distance modulus of m — M = 26.5, extinction of Af606W = 0-l ; 
and metallicity of [M/H] = —0.4, matching what we derived for the region surrounding NGC 300 
OT2008-1. We varied the IMF from a = -2.0 to a = -2.7 for a binary fraction of 0.35, and 
varied the binary fraction from 0.2 to 0.5 with a Salpeter IMF of a = —2.35. We ran MATCH on 
each of these models assuming a Salpeter IMF and binary fraction of 0.35, as in our original SFH 
derivation, to see if the age of the recovered burst changed with the input IMF or binary fraction. 
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We found that the only changes were in the amplitude of the burst, as expected since changing the 
IMF changes the normalization of the SFH, and that in each case MATCH still recovered a burst of 
8-13 Myr. Hence, even if the IMF or binary fraction is different from our assumed value, it should 
not affect our age or mass estimates. 



Discussion 



4.1. Implications for SN 2008S/NGC 300 OT2008-l-type Transients 

Both NGC 300 OT2008-1 and the similar transient SN 2008S have very highly obscured progen- 
itors that are bright point sources at mid-infrared wavelengths but are undetectable in the optical. 
The spectrum of SN 2008S is similar to a SN Iln spectrum, but many scenarios, including the 



outb urst of a star, are consistent with the low velocities implied by the narrow lines (jDessart et al 
20091 ). To dat e, two possible explana tions have been put forward to explain this class of lumi- 



nous transient. IThompson et al.l (|2008l ) suggest that the progenitors were heavily dust enshrouded 
AGB stars of mass ~ 10 Mq with a photospheric radius in the mid-IR of a few hundred AU, and 
that the transients may have been electron-capture SNe. iBotticella et al.1 (|2009l ) also endorses the 
electron-capt ure SN explanation for these transients, with a mass estimate of 6-8 Mq for SN 2008S. 
Alternatively, ISmith et al.1 (120091 ) argue that the transient may be due to a super-Eddington wind 
being driven from a > 20 Mq star. Th e SEP and light curve of M85 OT2006-1 are also very similar 
to NGC 300 OT2008-1 an d SN 2008S (|Kulkarni et al.ll2007l ). with a mass estimate on the low end, 
< 7 M Q (bfek et al.ll2~008h . 



The 12-17 Mq main-sequence turnoff mass we have derived as a lower limit for the precursor is 
in between the mass ranges suggested for bot h the electron-capture SN in an AGB sta r (6-11 



2008 



Bottice 



Smith et al 



la et al.l 120091 ) and the wind mechanism suggested by ISmith et al 
( 20091 ) suggests a rough lower limit of ~ 15 solar masses for the 



Tho mpson et al . „ 
(2009J). Although 

super-Eddington wind, there are no firm theoretical arguments that would exclude lower masses. 
Given both the uncertainty in these mass ranges and the possibility that our data is consistent with 
a mass as high as 25 Mq or as low as 11 Mq, we cannot at present rule out either option. 



Our mass estimate is consistent with the ~ 10-15 Mq estimate of iBond et al.l (J2009J) from 
Spitzer mid- to far-infrared fluxes. They infer an outflow from the spectrum of the t ransient, and 



propo se an OH/IR star which erupted and cleared its surrounding dust envelope. iBerger et al 



(J2009J) find evid ence of infall as we ll as outflow in their high-resolution spectra. Using a similar 
interpretation to ISmith et al.l (|2009l ). they suggest a blue supergiant or Wolf-Rayet star and give a 
more generous mass estimate of ~ 10-20 Mq, also fully consistent with our derived SFH. 
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4.2. The Red Supergiant Problem 



The progenitor mass for NGC 300 OT2008-1 reported in this paper has interesting implica- 
tions for stellar evolution theory. Prior to SN 2008S and NGC 300 OT2008-1, episodic mass loss 
evident in LBV outbursts and SN Hn was associated only with the most massive (> 20 Mq) stars 



(Smithetalj2004 



Gal- Yam et al 



capture SN (IPrieto et al.l l2008al: iT hompso n et a 



2007 



LBV eruptions (jSmith et al.l 120091 ; iBerger et al 



Smith 2007). Whether NGC 300 OT2008-1 is an electron- 



2008) or a non-disruptive outburst analogous to 
20091 ) . the narrow emission lines are a signature 



of significant mass loss at speeds lower than a typical SN. In either case, the association of a pro- 
genitor mass as low as ~15 Mq with an outburst and implied episode of mass loss is a new and 
somewhat unexpected result. Perhaps t he only obvious analog is SN 1987A, which was a ~ 2 



19881 : IWooslev 



1988 



Mq blue supergiant j ust before ex plosion (jWooslev et al.lll987l ; IWalborn et al.lll987l ; IWooslev et al 



Arnett Il99ll ). Its progenitor was likely a red supergiant (RSG) ~ 30, 000 yr 



earlie r ( Fransson et al. 19891 ). and may have made the transition via a LBV-like eruption (Smith 
20071 ). Might NGC 300 OT2008-1 be the signature of a similar transition? 



Whether significant mass-loss events are common for stars < 20 Mq has important conse- 
quences for our understanding of stellar evolution, especially during the last stages. If the cumu- 
lative effect of events like NGC 300 OT2008-1 is to expel a significant fraction of the hydrogen 
envelope, then we know that not all ~ 8-20 Mq stars experience such events, since the most com - 
mon type of SN, SN II-P, originate from stars in this mass range. However, ISmartt et al.l ((2009) 
have identified a potential lack of SN II-P progenitors more massive than 16.5 ± 1.5 Mq. RSGs, the 
progenitors of SNe II-P, should be common up to a maximum ma ss of at least ~ 25 Mq, since more 
massive stars will become Wolf-Rayet stars (e.g.. ICrowtheii 120071 ). Although statistic ally some SN 



proge nitors of masses 17-25 Mq should have been in the volume-limited sample of ISmartt et al 



(J2009J), they re port with 2.4a confid ence that RSGs in this mass range are not exploding as SNe II- 
P as expected. ISmartt et al.1 (J2009j) consider several resolutions to this problem, including the idea 
that stars in this mass range explode as another type of SN, but they favor the idea that the 
minimum progenitor mass for black hole formation is lower than expected, resulting in failed or 
weak SNe for this mass range. 

Intriguingly, the progenitor mass we find for NGC 300 OT2008-1 is consistent with the mass 
range that defines the onset of the SN II-P deficit, and the event itself could signify the transition 
away from a RSG phase. We consider whether stars in the mass range ~ 17-25 Mq experience 
optical transients (OTs) and episodic mass-loss, expelling much or some of their hydrogen envelopes 
and aborting a connection to SN II-P. Using the observed rate of SNe II-P and assuming the Salpeter 
IMF, we can estimate the minimum number of NGC 300 OT2008-l-like transients necessary to 
explain the RSG problem. If we assume that a minimum of one OT per star is necessary to keep 
it from exploding as a SN II-P, integrate the IMF from 17 to 25 Mq, and compare this to the IMF 
integrated from 8 to 17 Mq, this gives us a minimum OT-to-SN II-P ratio of 0.23 to explain the 
lack of high mass SN II-Ps. 
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Smartt et al.l (j2009l ) report 54 SN II-P in their volume limited (< 28 Mpc) survey during a 



span of 10.5 years. The deficit of higher mass SN II-P progenitors suggests that (in steady state) 
at least ~ 12 additional stars within the same volume must have shed some of their hydrogen 
envelopes in an OT during the duration of the survey; an even higher rate of OTs may be expected 
if multiple mass loss events are needed to fully strip the hydrogen env elope. Three simila r transients 
were reported between 2006 and 2008, one in 2006 and two in 2008 (jBerger et al.1 12003 ) . They are 
NGC 300 OT2008-1 at ~ 2 Mpc, SN 2008S at ~ 4 Mpc, and M85 OT2006-1 at ~ 1 7 Mpc. If we 



take t his as a rough rate of OTs, then the number of OTs over the > 10 yr span of the lSmartt et al 
(|2009l ) survey is ~ 11, and within a smaller survey volume. Therefore, the observed rate of OTs is 
not inconsistent with this type of event being sufficient to explain the paucity of high mass SN II-P 
progenitors b e tween 17 and 25 M@. However, given the uncertainties and small sample size in the 



Smartt et al.l (|2009i ) mass estimates, the different effective survey volumes for SN II-Ps and OTs, 



and the lack of a secure progenitor mass for SN 2008S, we cannot yet draw a firm conclusion. 



4.3. Limitations of the Method 

The work above demonstrates that one can obtain accurate precursor masses even when the 
actual precursor is undetected. Although the images we used were taken several years before the 
transient, there was nothing in the analysis that required direct identification of the precursor. 
Instead, the mass was derived through "guilt by association" using the properties of the underlying 
burst of star formation. In principle, one could use this method to derive masses of precursors of 
every transient or SN in the historical record. In practice, however, this method is somewhat more 
limited. 



4-3.1. Crowding 

Recovery of recent SFHs requires CMDs of individual main-sequence and helium-burning stars. 
However, images of individual stars are blurred together, with the crowding being more severe 
at larger distances, in high surface brightness regions, and/or for fainter stars (which are more 
numerous). Even with HST's resolution, crowding quickly limits the ability to measure fluxes for 
all but the brightest stars beyond ~ 4-5 Mpc. The crowding limit can be reached many magnitudes 
brighter than one would expect for a given telescope aperture and exposure time. This limitation 
is less of a problem for stars on the blue and red helium burning sequences, which are far more 
luminous than main-sequence stars of comparable masses. However, these sequences are only well- 
populated after ~ 25 Myr, making them useful only for low mass precursors (M < 9 Mq). Thus, 
crowding is a significant limitation on the applicability of this method, at least until the launch of 
telescopes that can surpass HST's angular resolution in the optical. Advances in infrared space 
telescopes may not bring much progress, given the faintness of main-sequence stars in the infrared. 
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4-3.2. Depth 

Even when crowding is not significant for a main sequence star of a given stellar mass, the SFH 
method requires depths that are a magnitude or more below the main-sequence turnofF of the stellar 
masses of interest. The recent star formation is largely constrained by the luminosity function of 
the main sequence, so that the behavior on the more sparsely populated upper main sequence is 
forced to be consistent with the numbers of stars on the lower main sequence. Thus, bursts are 
best constrained when the main sequence is well-measured substantially below the turnofF. 



4-3.3. Very Massive Stars 



While the modeling in this paper yielded excellent limits on the mass of this particular pre- 
cursor, we anticipate larger uncertainties for higher mass precursors, for several reasons. First, the 
optical luminosity of a very massive main-sequence star is a tiny fraction of its bolometric luminos- 
ity, making optical CMDs poor constraints on the masses of the most massive stars. The masses 
of the most m assive main sequence stars can only be accurately determined through spectroscopy 
(jMasseyl 120031 ) . given that stars with T > 30,000K have essentially identical optical and near-UV 
colors. Thus, while using SFHs to measure precursor masses may work for up to 20-30 Mq, it is 
less likely to be useful for higher masses. 

Second, the very youngest, most massive stars are likely to remain shrouded within their natal 
molecular clouds to some extent. If these stars have not emerged from their dust cocoons by ~ 5 
Myr (the turnofF age of a ~ 25 Mq star), then their measured luminosities and colors are likely to 
be significantly affected by dust, causing the main sequence to broaden to the red and to shiFt to 
lower luminosities. SFH recovery codes typically take dust into account statistically by including a 
"differential reddening" term to help match the width and position oF the main sequence. However, 
this correction cannot compensate For stars that are simply absent due to very high dust extinction; 
observations in the Magellanic clouds suggest that this eFFect is significant For ages less than 2 Gyr 
(jMassev et al.lll995l ). Dust can also be corrected For on a star-by-star ba sis, provided that imaging 
on the UV side oF the Balmer break is included (jRomaniello et al.ll2002l ). 



Third, very massive stars are likely to be Found within dense stellar clusters that have not yet 
experienced significant cluster dissolution. The timescale For this dissolution is ~ 10 Myr , such 



that 50% oF stars w ith mass 15 Mq are likely to still be bound into clusters (jLada Lada 



2003 



Bastian et al.ll2005l ). These clusters have higher than average crowding, and thus stars within them 
are prone to larger photometric errors and brighter limiting magnitudes. 

Finally, the upper main sequence is always less populated than the lower, making it statistically 
possible that the youngest apparent "turnofF" in the main sequence appears to be older than the 
true turnoff, due to the sparse sampling oF the high mass end oF the IMF. 
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4-3.4- Region Size 

An additional concern with the method is the choice of region size. If the radius within which 
stars are selected is too small, it is possible that stars will have dispersed from an unbound cluster 
to beyond this radius. However, the dispersal timescale is tens of My r, as the stars are expecte d 



to disperse with the typical cluster velocity dispersion of a few km s _1 ( Bastian fe Goodwinll2006l ). 
In the 5" (50 pc) region we selected, if the cluster velocity dispersion is 1 km s -1 = 1 pc Myr -1 
(a reasonable value for a small cluster), stars would still be contained within the region after 50 
Myr, which is longer than the estimated timescale of star formation for this cluster. For a larger 
cluster with a velocity dispersion of 3 pc Myr -1 , the stars would remain within 50 pc for 17 Myr, 
still longer ago than the beginning of the star formation event we see in this region. 

Choosing a larger radius would ensure that all stars from the cluster are included; however, 
a larger radius also increases contamination from surrounding unrelated stars. We find that using 
a radius of 7.5" instead of 5 (2.25 times the area enclosed) results in a less well-defined burst, 
probably because we are enclosing stars slightly older than the population which gave rise to the 
transient. 



4-3.5. Choice of Stellar Evolution Models 



We have used the isochrones from the Padova group (jMarigo et al.l l2008l ) in our analysis, 
and the errors reported assume that those models are correct. How e ver, a dditional systematic 
errors may be introduced by uncertainties in the models. iGallart et al.l (|2005l ) provide an in-depth 
discussion of the discrepancies between different models and their effect on deriving SFHs using 
CMD fitting, and they conclude that the uncertainties seem to be smallest for young stars, which 
is reassuring for studies of transient progenitors. To assess the e rrors in more detail, we c ompare 



the Padova models with the currently available Geneva models (ILejeune fe Schaererll200ll ). Both 



sets of models include the effects of convective core overshooting, but not rotation. The Padova 
models include thermal pulses on the AGB, while the Geneva models do not. 

Figure [6] shows isochrones from the two models for ages 8 Myr and 13 Myr, the flanking ages of 
the star formation event we found for NGC 300 OT2008-1. Metallicity is Z = 0.008 ([M/H] = -0.4). 
The left panel plots the isochrones for both models for each age, with the main-sequence turnoff 
magnitude marked. The models are very similar up to the turnoff, and diverge for more massive 
stars. We are especially concerned with the mapping between age and stellar mass, since we use 
this relationship to go from an age estimate to a stellar mass estimate for the NGC 300 OT2008-1 
progenitor, so we plot mass as a function of absolute magnitude for the same isochrones in the right 
panel of Figure [6l For each age, the mass-luminosity relationship is identical between the models 
up until the main-sequence turnoff, but there are discrepancies above the turnoff. Since the age of a 
young stellar population is primarily determined by the main sequence, we expect that differences 
in the models related to the treatment of post-main-sequence phases will not significantly affect 
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our results. 



A potentially significant effect not addressed in the previous comparison is rotation. For solar 
metallicity and an initial rotation velocity of 200 km s -1 , the age esti mate for a star of a give n 
luminosity would increase by about 25% from the non-rotating model (jMeynet &: Maederl 12000 ) . 
Hence, including rotation in the models would increase our age estimate for the star formation 
event to 10-16 Myr and the upper mass limit for the turnoff would decrease to 15 Mq. The lower 
mass limit for the turnoff changes from 12.3 Mq to 11.9 Mq. Hence, even including the effects 
of rotation at this velocity, our turnoff mass range of 12-17 Mq is still vali d. Rotation sp e eds in 
excess of 200 km s~ 1 could lead to smaller possible masses for the progenitor. IVazquez et al.1 ( 20071 ) 
show that the decrease in main sequence lifetime also holds for stars at Z = 0.008 (the approximate 
metallicity of the region around the transient), but isochrones for this metallicity are only available 
for masses > 30 Mq. 



4-3.6. Ambiguous Bursts 

The final possible limitation is that Nature may not always be as kind as she has been in 
this particular instance. While NGC 300 OT2008-1 came from an unambiguous isolated burst, 
some transient phenomena may be associated with regions containing multiple bursts at a range 
of ages, making the association of mass ambiguous. This effect is more likely to be pronounced 
in more distant galaxies, where a typical aperture size of a few arcseconds encompasses a larger 
physical area of the disk, and thus is more likely to contain more than one dissolving open cluster. 
On the other hand, t he analysis of star formation regions in M81 (D ~ 3.7 Mpc) presented in 



Gogarten et al.l (120091 ) uses ~ 27" « 500 pc regions, and finds bursts that are well localized in time 



up to 100 Myr. 



4.4. Benefits of the Method 

In spite of the above limitations, there are a number of clear benefits to the procedure we 
have demonstrated in this paper. The first and most obvious one is the removal of the need for 
precursor imaging in determining the mass. Instead, targeted observations after the event are 
sufficient, and can be optimized for recovery of the SFH, through proper choices of filters and 
observing strategy. Second, unlike in precursor imaging, one does not need accurate astrometry. 
Since the SFH recovery uses stars within a few arcsecond wide aperture, it is not required that the 
SN or transient be localized to better than half the width of the selection aperture. Third, given the 
uncertainties in the late stages of stellar evolution, it is perhaps easier to derive an accurate mass 
from the bulk properties of the better-understood main- and helium-burning sequences, rather than 
from the luminosity and SED of a star on the brink of explosion. Finally, the relationship between 
progenitor mass and SN type as predicted by stellar evolution models relies on assumptions about 
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mass loss during late stage stellar evolution. Contributions to the existing catalog of precursor 
masses using the SFH technique will improve the statistics in comparing observations with stellar 
evolution theory, and may even give more reliable initial mass estimates than can be determined 
from an individual star. 



5. Conclusions 

In this paper we have demonstrated a technique for measuring the masses of precursors of 
luminous transient events, even when no pre-event imaging exists. By using standard tools to 
recover the SFH from stellar populations, we find that NGC 300 OT2008-1 originated from stars 
formed in a burst between 8 and 13 Myr ago with a certainty of 70%. The current turnoff mass 
associated with this burst is 12-17 MQAssuming the transient is due to an evolving massive star, 
then the mass of the precursor must be higher than this turnoff mass, but less than the 16-25 Mq 
mass limit above which stars of this age have exploded. The res ulting mass range of 12-25 M P) 



agrees quite well with estimates of 10-15 M bv lBond et al.1 (|2009l ) and 10-20 M bv lBerger et al 



(|2009l ). This technique therefore shows great promise for significantly expanding the number of 



SNe and transients with reliable precursor masses. 

The SFH method presented here for determining precursor masses is largely complementary 
to direct precursor imaging, as each method has benefits and limitations. At present, the limits 
of the SFH technique are largely observational. The angular resolution of HST currently restricts 
this technique to well within D < 10 Mpc. In a future paper (Gilbert et al. 2009, in preparation), 
we will present an application of the method to all SNe and transients for which sufficiently deep 
data exists in the HST archive. 

With future space-based optical telescopes one can certainly push the technique to much larger 
distances, greatly expanding the volume of candidate SNe. The more intensive monitoring of nearby 
galaxies can also significantly increase the number of S Ne for which t his technique can be used; 



thankfully such monitoring is currently underway (e.g., iLeamanl 120081 : iDi Carlo et al.ll2008l ). The 
main benefits of the method are that imaging can be done after a transient event, even without 
accurate astrometry, and that fitting the entire main sequence of a star formation event may provide 
more reliable mass estimates. 

In closing, we note that precursor imaging is still highly desirable. The technique we have 
employed gives only a constraint on the main-sequence mass of the star that eventually erupted. 
However, it gives no constraint on the exact phase of stellar evolution that the star was in im- 
mediately before the eruption. Precursor imaging can yield much more information about a star 
than its mass alone (e.g., color, magnitude). Given the impact of the unexpected phase of the 
SN 1987A precursor (i.e., that it was a blue, rather than red, supergiant at the time of explosion) 
on our understanding of core-collapse SNe, any additional information on the precursor is likely to 
be highly significant. 
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Fig. 2. — Left: recent star formation history (< 100 Myr) of the selected region in NGC 300, 
showing a significant star formation event at 8-13 Myr. SFH was found by comparing the observed 
ANGST CMD with synthetic CMDs based on theoretical isochrones. Error bars are the quadrature 
sum of the uncertainties from distance and extinction and the 68% confidence int erval from Monte 
Carlo simulations. Right: main-sequence turnoff mass vs. age (jMarigo et al.ll2008l ). The boundaries 
of the star formation event and corresponding stellar masses are marked. These masses represent 
the lower limits on the mass of the transient. 
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Fig. 3. — Simulated bursts of age 8-13 Myr and varying amplitudes (gray shaded regions) and 
derived SFHs. The cumulative distribution functions show the total star formation summed from 
older ages to younger ages. The central black line is the median value of the Monte Carlo tests, and 
the surrounding hashed region is the la error range. Derived SFH and errors for NGC 300 OT2008- 
1 are in red. The number of stars refers to upper main-sequence stars (F606W — F814W < 0.3, 
F606W < 26). 
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Fig. 4. — Simulated bursts with ~ 50 upper main-sequence stars and varying durations (gray 
shaded regions) and derived SFHs. The central black line is the median value of the Monte Carlo 
tests, and the surrounding hashed region is the lex error range. Upper row: constant end time 
(bursts start at 10, 16, and 20 Myr ago and end at 8 Myr); lower row: constant start time (bursts 
start at 13 Myr ago and end at 5, 6, and 10 Myr). Burst duration is labeled on each plot. 
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Fig. 5. — Simulated bursts with ~ 50 upper main-sequence stars and varying ages (gray shaded 
regions) and derived SFHs. The central black line is the median value of the Monte Carlo tests, 
and the surrounding hashed region is the la error range. Upper row: burst ages are 4-8, 5-10, 
6-10 Myr. Lower row: burst ages are 10-16, 16-20, and 20-25 Myr. Burst duration is labeled on 
each plot. 
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6. — Compariso n of the Padova ([Marigo et al.l 120081 . solid lines) and Geneva 



(jLeieune fc Schaererll200ll . dotted lines) isochrones for Z = 0.008 and ages of 8 Myr (blue) and 13 
Myr (red). The left panel shows the isochrones on a CMD, while the right panel shows stellar mass 
as a function of absolute magnitude. The magnitude of the main-sequence turnoff is highlighted 
with a dashed line for each age. The maximum masses shown by the endpoints of the isochrones 
cover a range of 16-25 M®. The two different isochrone sets agree well for main sequence stars, 
but show significant differences in their post-main-sequence evolution. 



